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Abstract 

Linear clusters made by tightly connecting two or more metallic nanoparticles have new 
types of surface plasmon resonances as compared to isolated nanoparticles. These new reso- 
nances are related to the size of the junction and to the number of interconnected particles and 
have direct interpretation as eigenmodes of a Boundary Integral Equation (BIE). This formula- 
tion allows effective separation of geometric and shape contribution from electric properties of 
the constituents. Results for particles covered by a thin shell are also provided. In addition, the 
present analysis sheds a new light on the interpretation of recent experiments from literature. 

The interaction of electromagnetic fields with discrete systems ranging from plasmonic nanopar- 
ticles 1,2 to biological cells^ is governed by their polarizabilities and is strongly dependent on both 
dielectric and geometric properties of particles and surrounding medium. Metallic nanoparticles 
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such as noble-metal nanostructures exhibit the localized surface plasmon resonance (LSPR) phe- 
nomenon when the electromagnetic radiation interacts with the collective oscillations of the con- 
duction electrons. 1 The resonances in visible (VIS) and near-infrarred (NIR) induce extremely 
strong fields confined to the proximity of the surface of metallic nanoparticles, and therefore cre- 
ating the opportunity for highly sensitive sensors, 4 surface enhanced Raman scattering (SERS), 5 
surface enhanced infrared scattering (SEIRS), 6 near-field microscopy, 7 photoluminescence,- sec- 
ond or higher harmonic generation,-^^ and other applications. Closely spaced or nearly touching 
(connected) nanoparticles— - — as well as clustered particles in different configurations 14 provide 
additional enhancements of the field within the junction region. Moreover, the "dumbbell" shaped 
particles have also shown other remarkable properties like the large red-shift in the dipolar re- 
sponse,-^^ and large variations in the spectral response.-^*^ 

We demonstrate the spectral BIE approach by calculating the dielectric response for clusters 
made of up to four rotationally symmetric particles of two kinds: elongated particles, longer along 
the rotation axis (aspect ratio greater than unity) and flatten particles, shorter along the rotation axis 



(aspect ratio less than unity). [figure][l][] 1 shows the imaginary part of the polarizability for such 
clusters, when the electric field is along the rotation axis. Resonances in the mid-infrared (mid-IR) 
region of the spectrum appear for both kinds of clusters and are associated with the joints between 
the particles. These resonances are enhanced as the relative size of the joints become smaller and 
they. 3 The additional resonance for clusters of more than one flatten particles (sub-figure b) appears 
because the cluster as a whole becomes elongated. This resonance is visible for clusters made of 
any number of elongated particles (sub-figure a). There are no new resonances introduced when 
the number of particles increases. In recent works-^^*^ the distinct optical behavior of touching 
dimers has been attributed to the joints as soon as the particles come in conductive contact in 
the dimer. 13 Our results indicate that the resonances due to the joints are located far in mid-IR 
region of the optical spectrum. On the other hand, in the papers cited above,- 1 ^^ the shape of 
particles used in numerical simulations was close to spherical shape, which is the border shape of 



flatten particles. This observation together with the results of [figure][l][]l suggests that the red- 
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shifted plasmon resonance observed in-^i^ for dimers is the additional resonance which appears 
in clusters made of flatten particles. 

The idea of treating surface plasmons as electrostatic resonances has been introduced in— by 
showing that source-free electromagnetic fields correspond to the eigenvalues of the BIE. 17 Using 
a similar approach, we show that a spectral representation of BIE provide the dielectric responses 
of homogeneous nanoparticles and nanoshells in compact forms, which can be directly used in 
designing plasmonic structures. The design of analyte sensing is given as an example, where the 
explicit equations of plasmon resonances are used to calculate the refractive index sensitivity for 
each eigenmode. 
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Figure 1 : Imaginary part of polarizabily for clusters made of elongated (a) and flatten (b) particles 
in longitudinal fields. Mid-IR resonances are associated with the interparticle joints and become 
more prominent as the joints tighten. Resonances in the visible spectrum are associated with the 
shape and with the number of particles. 
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Each eigenmode of the linear response operator has a contribution in the expression of polar- 
izability proportional with its coupling weight and inversely proportional with the depolarization 
factor. Although this property has been demonstrated for biological cells, 3 it is also valid for metal- 
lic nano-particles and nanoshells. If an eigenvalue has a non-zero weight and is close to the value 
1/2, then it is responsible for a mid-infrared (mid-IR) plasmon resonances, which can be used in 
optical nanoantennas.— ^ As the interparticle junction tightens, the eigenvalue approaches 1/2 and 
the resonance is enhanced. The coupled-particle LSPR occurs at a frequency red-shifted from that 
of a single-particles. 2 

The nanoparticle clusters are considered under the electrostatic (quasistatic) approximation,— 
when their size is no more than one tenth of the wavelength of the incoming radiation. The op- 
tical behavior of metallic nanoparticles 16,17 as well as the dielectric behavior of living cells 3 can 
therefore be obtained from a solution of the Laplace equation derived from the corresponding BIE 
written for the surface E which defines the particle. The polarizability of such an object of permit- 
tivity £i embedded into a medium of permittivity £ is given in,- 3 -^ 7 - and has the following spectral 
representation 

Here electric parameter A = (£; — £ ) /(£,• + £ ), while %k an d Pk are eigenvalues and coupling 
weights corresponding to the k-th eigenmode of the linear response operator M. This geometric 
operator depends only on the shape of the particle, is defined by 

and have the following properties:- 1 ^ its spectrum is discrete, real and bounded to the interval [- 
1/2, 1/2]; the value 1/2 is always an eigenvalue, and it has only one active eigenvalue for spherical 
particles. The imaginary part of the polarizability a is directly related to optical absorption of 
metallic nanoparticles,— and depends on the polarization of the incident radiation 3,17 through 
dipole-coupling weights p^. Equation 1 shows how material properties, all collected within A, are 
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separated from geometric properties contained in %k an d Pk- 

An object covered by a very thin shell of permittivity £5 has a polarizability calculated from ?? 
by replacing X by X k = (£ k — E )/{E k + £ ), 3 - where the effective mode-dependent permittivity £& 
is given by: 

/ £i~£s \ 

£k = £s { 1 + e s + 8{\/2- Zk )£ i + 8(1/2 + Xk )e s ) (3) 
where 8 is a small parameter roughly equal to the relative volume of the shell with respect to the 
particle. 3 This equation allows a direct connection to the hybridization model 18 

Within various ranges of frequencies, the dielectric functions have in general the form £ = A + 
B/W, where A and B are nearly constant functions of angular frequency 00, while W is a fast varying 
function of CO. In the radiofrequency domain £ = £ — io/(e vac co) with £ the static permittivity, 
£ vac the vacuum permittivity and o the DC conductivity. In the optical regime dielectrics can 
be successfully described by a nearly constant dielectric function £ = Ed and metals by a Drude 
dielectric function £ = £ m — 60%/ G)/((Q + iy), where e m is the interband contribution, cop is the 
electron plasma frequency and 7 is the damping constant. 

The frequency-dependent polarizability of a homogeneous metal nanoparticle of arbitrary shape 
can then be written as a sum of Drude-Lorentz terms by using expansion ??: 

™ f,^ V Pk( e m — E-do) rA , 
ttplasmon{CO) = > " (4) 

k £ eff 
Pk £do ®pk 



1/2 - Zk e e ff (o{(o + iy) - (b 2 pk 



where E& is the permittivity of the outer dielectric medium and the effective dielectric parameter is 
£ e ff = ( 1/2 + %k)Edo + ( 1/2 — %k)E m - For very small particles/clusters dissipation can be neglected 
(7 <C (Op) and G) p k is the frequency of the localized plasmon resonance, which depends explicitly 
on the dielectric and geometric parameters of the particle: (b 2 pk = ( 1 /2 — %k) °tf> I eef f ■ 

Each active mode in ?? has a fast-varying (Drude-Lorentz term) factor proportional to its cou- 
pling weight p k and inversely proportional to (1/2 — Xk)> me depolarization factor. In clustered 
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particles with tight junctions, several eigenvalues get closer to 1/2 and the factor pk/ (1/2 — Xk) 
has an important contribution even though the couplings pk are vanishingly small. The local field 
enhancement factor (the height of the resonance peak) is Pk/(l/2 — Xk)(£do/ £ e ff ) ( ®pk I Y) as m tne 
harmonic oscillator model.— As Xk approaches 1/2 for tighter and tighter junctions the plasmon 
resonance red- shifts toward lower frequencies. 

The polarizability of dielectric particles of permittivity covered by a thin metallic layer and 
embedded into a dielectric medium of permittivity e^ is obtained in the first order in 8, from ?? 
and ?? as 

"shell — (5) 

k fc eff 

Pk £do / W pk 



1/2- Xk e e ff V o)( co + if) - co' 1 



pk 



w P k 



cd(co + iy)-co'V-co 2 /£ n 

with fig = (5(1/2 + Xk)(l/l-Xk)a$)/e ftS , co'pl = (e di /e m ){8{l/2 -Zfc) 2 ^)/e e ff and e eff = 
(1/2 + Xk) £ do + (1/2 — Xk) £ di- Like me hybridization model, 18 ?? exhibits two plasmon reso- 
nances. One resonance frequency for a shelled particle is red-shifted with respect to the one of the 
homogenous particle by a factor of ^8(1/2 + Xk)- Th e value of the peak of the resonance is also 



reduced by the factor \/8{\/2 + Xk)- Results obtained in 19 can be derived from ??, by observing 
that £ e ff determines the influence of e^i and e^ on the plasmon resonance in nanoshells. The other 



resonance, located at ^ co'^ + CO 2 / e m is blue-shifted with respect to co' pk and (o p k and is smaller 
than both the red-shifted resonance and the resonance of a similar homogeneous particle. 

A straightforward application of ?? and ?? is in the design of analyte sensing. Large changes 
in the position of the plasmon resonance in response to small changes in the refractive index of 
the outer medium are desired. The refractive index rido sensitivity for the wavelength k p k of each 
dipole active eigenmode is dX p kjdn do = (1/2 + Xk) n do^pk/ e eff which increases with increasing 
wavelenght and Xk- A body of theoretical and experimental research 2 has suggested that either 
elongated particles or thin nanoshells have high refractive index sensitivity. A direct explanation 



of this derives directly from ?? and ?? which show the red shift of the plasmon resonance for 
tightly joined clusters, elongated particles and nano shells. A figure of merit is given by the 
ratio of refractive index sensitivity over the width of the plasmon resonance. In the limit 7 <C CO p , 
the resonance width is given by the damping constant 7. In the case of nanoparticles 7 is given 
by the the damping constant of the bulk and an additional term due to the surface scattering of 
conduction electrons: y = Vf(l/l + l/L),— with Vf the Fermi velocity (about 1.4 x 10 6 m/s for 
both Ag and Au), / the bulk room temperature mean free path (about 52 nm for Ag and 40 nm for 
Au), and L is the effective mean free path for collision with the particle's boundary. For a spherical 
particle of diameter d, L = 2d/3, such that if d = 25 ~ 30 nm, 7 is about 0.1 eV. The damping 
constants in nanoshells are even greater due to larger effective mean free paths for collision;— 
therefore nanorods are better than nanoshells when used in analyte sensing. Using ?? and ?? the 
refractive index sensitivity of nanoparticles can be directly calculated as a function of the aspect 
ratio or as a function of composition of the nanoparticle clusters. 

We have analyzed clusters of nanoparticles as n touching objects with a rotational surface 
generated by the equation: x = g(z) cos (j>,y = g(z) sin0 and — L n <z<L n (the insets of Figure la 
and lb) with 

g(n,h,A,a,b,z) = A(l + sign(z sign(z) - (n-l)a)) x 
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1 — ((z — sign(z) (n—\)a)/a) 
1 + (1 + b(z- sign(z)(n - \)a) 2 



1 — floor 



/ sign(z) 
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l + (l-#(z)) 2 



where H(z) = mod(( — i)fl° or (z/ a )+' 1 - 1 z,a) 2 , sign(z) is signum function, -1, or 1 if z is negative, 
zero or positive, floor(z) is the greatest integer less than z, and mod(x,y) is the remainder of the 
division of x by y. Parameters A, a, b determine the shape of an individual particle, for example, 
the ratio A/ a is the aspect ratio, and h determines the size of the connecting gap, and L n = na. 



[figure] [1][]1 shows the resonances in the mid-IR introduced by tight junctions because some 



eigenmodes have eigenvalues near 1/2 and non-zero, but very small, weights. In a recent paper we 
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Figure 2: Imaginary part of polarizabily for clusters of elongated (a) (A = 2, a = 3, h = 0.2, 
b = 0. 1) and flatten (b) (A = 3, a = 2, h = 0.05, b = 0. 1) particles in visible spectrum. 
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have shown that in a cluster of n connected particles there are n — 1 eigenvalues close to 1/2 (the 
smaller the junction the closer the eigenvalues). 3 However, only the eigenmodes with antisym- 
metric combination of charge distribution in the lobes are active and contribute to polarizability. 
Thus the eigenmodes close to 1/2 induce mid-IR plasmon resonances in optical spectra of metal- 
lic clusters even though the weight (dipole coupling) is quite small. In our examples clusters of 
two and three particles have just one resonance, while for clusters of four particles there are two 
overlapping resonances. As a rule, all resonances move toward infrared with increasing number of 



particles in the clusters, as seen in [figure] [2][]2 
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Figure 3 : Imaginary part of polarizabily for clusters of elongated (a) and flatten (b) particles placed 
in electric perpendicular to the rotation axis. 



In contrast to parallel fields ( [figure] [1][]1 and [figure] [2][]2| ), the polarizability for clusters in 



field perpendicular to their axes, displays no new significant resonances, as seen in [figure] [3] []3 
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The existing resonances move toward ultraviolet as n increases. Also, the mid-IR resonance shifts 
toward larger wavelengths as the ratio h/A decreases. 



[figure] [2] [] 2 and [figure] [3] []3 indicate that, in clusters as well as in individual nanoparticles, 



the amplitude of plasmon resonance depends also on the direction of polarization. Thus flatten 
particles have the greatest plasmon resonances for transverse polarization, while elongated parti- 
cles has theirs for longitudinal fields. The explanation is based on the harmonic oscillator model, 
which relates the amplitude of plasmon resonance to the polarizability of the particles.— 

To conclude this letter, the spectral BIE is a rational approach which provides useful insights 
of the plasmonic behavior of metal nanoparticles. Linear clusters of overlapping nanoparticles 
are good candidates for optical sub-wavelength guides.- 1 ^ 3 - Our results suggest that the collective 
resonances originating from single-particle resonances can be easily activated by exciting just one 
particle at one end of the cluster. In contrast, the collective resonances that are not present in an 
individual particle can be successfully activated by exciting the last two particles at the end of the 
cluster. Another possible application of linear clusters of overlapping nanoparticles is their use as 
optical nano-antennas instead of nanorods.— "21 The utilization of linear clusters as optical nano- 
antennas has a lot of promise due to recent advances in the manufacturing of highly uniform dimer 
arrays with gaps smaller than 10 nm. 25 Thus instead of using long and narrow nanorods, one can 
use linear clusters, with lengths comparable to the length of the nanorods and junctions comparable 
to the cross section of the nanorods. 3 
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